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ABSTRACT 


One  hundred  seven  commercially  available,  off-the-shelf,  1 .2-Ah  cells  were  tested  for  calendar  life  and  CS 
cycle-  and  CD  cycle-life  using  the  new  USABC  PHEV  Battery  Test  Manual.  Here,  the  effects  of  temperature 
on  calendar  life,  on  CS  cycle  life,  and  on  CD  cycle  life;  the  effects  of  SOC  on  calendar  life  and  on  CS  cycle  life; 
and  the  effects  of  rest  time  on  CD  cycle  life  were  investigated.  The  results  indicated  that  the  test  procedures 
caused  performance  decline  in  the  cells  in  an  expected  manner,  calendar  <  CS  cycling  <  CD  cycling.  In  some 
cases,  the  kinetic  law  changed  with  test  type,  from  linear-with-time  to  about  ^.Additionally,  temperature 
PHEV  battery  testing  was  ^oun(^  t0  stress  the  cells  more  than  SOC,  causing  increased  changes  in  performance  with  increasing 

Performance  degradation  temperature. 
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1.  Introduction 

Lithium-ion  batteries  continue  to  attract  much  interest  in  appli¬ 
cations  where  high  specific  or  volumetric  power  or  energy  is 
required.  High-energy  lithium-ion  batteries  are  also  being  con¬ 
sidered  for  automotive  applications  by  the  U.S.  Department  of 
Energy-supported  U.S.  Advanced  Battery  Consortium  (USABC)  [1]. 
The  batteries  usually  consist  of  a  metal-oxide  positive  electrode, 
a  carbon  negative  electrode,  and  an  organic  electrolyte  containing 
dissolved  lithium  salts. 

Layered-oxide  cathodes,  such  as  Li(Ni1/3Mn1/3Co1/3)02,  repre¬ 
sent  good  candidates  for  automotive  applications  because  of  their 
high  specific  capacity  [2].  However,  the  cycle  life  of  this  oxide  was 
not  as  high  as  desired  [3].  On  the  other  hand,  spinel  oxides,  such  as 
LiMn204,  are  also  viable  candidates.  These  oxides  are  low-cost  and 
have  fast  kinetics,  which  makes  the  oxide  suitable  for  high-power 
applications  [2].  The  spinels  were  reported  to  have  lower  specific 
capacity  than  the  layered  oxides  [2]  and  to  degrade  rapidly  due  to 
manganese  dissolution  [3-6].  There  have  been  reports  in  the  litera¬ 
ture  that  blending  these  two  materials  produced  a  composite  with 
the  benefits  of  both  [7-12]. 
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Two  U.S.  Department  of  Energy  national  laboratories,  Argonne 
National  Laboratory  (ANL)  and  Idaho  National  Laboratory  (INL), 
continue  to  collaborate  to  understand  the  causes  of  performance 
decline  in  lithium-ion  batteries.  Results  from  this  collabora¬ 
tion  using  three  positive  electrode  materials— LiNio.8Coo.2O2, 
LiNio.8Coo.1Alo.1O2,  and  Lii.05(Mn1/3Co1/3Ni1/3)0.95O2-are  given  in 
Refs.  [13-18]. 

The  procedures  outlined  in  the  USABC  test  manuals  [19-21] 
are  intended  to  show  the  promise  of  a  technology  versus  a  set 
of  performance  and  cost  targets.  No  knowledge  of  the  actual  bat¬ 
tery  chemistry  is  needed.  Thus,  the  evaluation  of  cells  concentrates 
on  their  performance  and  life  and  how  their  life  is  affected  by 
SOC,  time,  temperature  and  type  of  test.  Calendar  and  cycle  life 
tests  were  performed  to  determine  the  aging  characteristics  of 
the  blended  cathode  material  under  plug-in  hybrid  electric  vehicle 
(PHEV)  testing  conditions  [19].  It  should  be  noted  that  the  calendar 
tests  in  the  PHEV  and  hybrid-electric  vehicle  (HEV)  [20]  manu¬ 
als  are  essentially  the  same.  The  cycle  life  test  is  more  complex, 
depending  on  what  is  to  be  learned.  Instead  of  just  one  cycling  mode 
of  operation,  charge-sustaining,  as  in  an  HEV,  the  PHEV  operates  in 
charge-depleting  (CD)  as  well  as  charge-sustaining  (CS)  modes,  as 
shown  schematically  in  Fig.  1 .  In  the  CD  mode  of  operation,  the  bat¬ 
tery  powers  the  vehicle  directly;  the  internal-combustion  engine  is 
not  used  at  all.  After  a  while,  the  battery  energy  becomes  exhausted. 
At  this  low  state  of  charge,  the  PHEV  will  operate  in  CS  mode,  similar 
to  that  of  a  hybrid-electric  vehicle. 


10214 


J.  Belt  et  al. /Journal  of  Power  Sources  196(2011)  10213-10221 


Time 


Fig.  1.  Schematic  showing  the  two  PHEV  operating  modes,  charge-depleting 
and  charge-sustaining  (after  Ref.  [19]).  In  charge-depleting  mode,  the  internal- 
combustion  engine  is  not  used;  energy  from  the  battery  alone  powers  the  vehicle. 
The  charge-sustaining  mode  is  similar  to  that  of  the  HEV;  the  power  from  the  bat¬ 
tery  supplements  that  from  the  engine.  The  abbreviation,  UECs,  shown  in  the  figure, 
is  the  SOC  window  from  which  the  usable  energy  in  CS  mode  is  based. 


The  PHEV-related  test  methods  and  protocols  are  new  and  have 
not  been  validated.  The  objective  of  work  below,  in  part,  is  to  val¬ 
idate  that  the  test  methods  stress  the  battery  as  expected.  In  the 
study  described  below,  the  effects  of  temperature  on  calendar  life, 
on  CS  cycle  life,  and  on  CD  cycle  life;  the  effects  of  SOC  on  calendar 
life  and  on  CS  cycle  life;  and  the  effects  of  rest  time  on  CD  cycle  life 
were  investigated  using  commercially-available,  18650-size  cells. 
The  results  of  this  work  will  help  elucidate  the  additional  stresses 
on  the  battery  when  it  is  cycled. 

2.  Experimental  method 


Fig.  2.  HPPC  profile  and  voltage  response  for  a  typical  cell. 
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2.1.  Testing 

One  hundred-seven  commercially-available,  off-the-shelf,  1.2- 
Ah,  18650-size  cells  were  used  in  this  work.  These  cells  contained 
a  physically-blended  Li-Ni-Mn-Co  layered-oxide  and  Li-Mn-0 
spinel  cathode,  an  organic  electrolyte,  and  a  carbon  anode.  The  test 
matrix  is  presented  in  Table  1.  All  cells  were  charged  using  the 
manufacturer’s  recommended  algorithm:  charge  at  1.2  A  to  4.2  V, 
followed  by  a  potentiostatic  hold  at  4.2  V  for  a  total  charge  time  of 
2  h  or  until  the  current  drops  below  50  mA. 

Before  the  start  of  the  aging  tests,  all  cells  were  characterized 
in  terms  of  their  10-kW  rate  capacities  and  by  the  hybrid-pulse 
power  characterization  (HPPC)  test  at  30  °C.  The  cell  groups  were 
aged  as  shown  in  Table  1  and  the  performance  was  averaged  for 
each  group.  The  cells  tested  at  temperatures  higher  than  30  °C  were 
heated  to  the  test  temperature  and  allowed  to  equilibrate  for  8  h 
before  the  aging  period  began.  Every  32  days,  testing  was  stopped 
and  the  cells  were  allowed  to  rest  at  30  °C  for  at  least  8  h.  Changes  in 
cell  performance  were  measured  by  repeating  the  characterization 
tests  at  30  °C  (a  reference  performance  test,  or  RPT).  Testing  was 
then  resumed  for  a  total  of  10  RPTs. 


Fig.  3.  Charge-sustaining  cycling  profile  [19].  The  power  levels  given  on  the  vertical 
axis  are  divided  by  the  BSF  to  yield  the  actual  power  levels  used  to  test  a  given 
battery.  The  unsealed  profile  removes  50  Wh  during  the  discharge  portion  of  the 
cycle;  the  regen  portion  charges  the  battery  with  46.2  Wh  of  energy,  producing  a 
net  removal  of  3.8  Wh. 


The  HPPC  test  was  a  constant-current  test  and  consisted  of 
removing  10%  of  the  rated  capacity  at  the  C/1  rate,  resting  for  1- 
h  and  applying  the  test  profile.  This  test  sequence  was  repeated 
for  a  total  of  9  profiles.  Fig.  2  shows  the  HPPC  test  profile  and  the 
voltage  response  of  a  typical  cell.  For  the  cells  tested,  the  discharge 
current  in  the  HPPC  profile  was  5.3  A,  and  the  regeneration  (regen) 
current,  4.0  A. 

Based  on  the  initial  HPPC  results,  the  average,  calculated 
battery-size  factor  (BSF)  was  1400  cells.  The  BSF  was  used  to  scale 
the  CS  and  CD  cycle  profiles  to  accommodate  the  capabilities  of 
the  test  cells.  The  generic  CS  and  CD  cycle  profiles  [19]  are  given 
in  Figs.  3  and  4,  respectively.  In  this  experiment,  all  SOCs  for  aging 
were  defined  in  terms  of  cell  potentials:  3.69,  3.89,  and  4.09  V  for 
30,  60,  and  90%  SOC,  respectively.  For  CS  cycling,  the  battery  was 
discharged  to  the  target  SOC  at  the  C/1  rate  and  allowed  to  rest  for 


Table  1 

Test  matrix  used  in  this  work.  Each  entry  in  the  table  represents  the  number  of  cells  tested  under  that  condition. 


Test 

Temperature,  °C 

SOC 

Rest  time,  min 

30 

40 

50 

60 

30 

90 

0  20 

40 

Calendar  life  (60%  SOC) 

Calendar  life  (30  °C) 

10 

5 

5 

4 

5 

5 

CS  cycle  life  (60%  SOC) 

CS  cycle  life  (30  °C) 

9 

5 

5 

5 

5 

5 

CD  cycle  life  (15-min  rest  time) 

CD  cycle  life  (30  °C) 

9 

5 

5 

5 

5  5 

5 
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Fig.  4.  Charge-depleting  cycling  profile  [19].  The  power  levels  given  on  the  vertical 
axis  are  divided  by  the  BSF  to  yield  the  actual  power  levels  used  to  test  a  given 
battery.  The  unsealed  profile  removes  545.61  Wh  of  energy;  about  6.2  repetitions  of 
the  profile  remove  3.4  kWh. 


1  h.  The  CS  profile  was  repeated  continuously  for  32  days  at  the  test 
temperature. 

Starting  from  the  initial  state  of  charge  (~90%  SOC,  the  PHEV 
charge  target),  the  CD  cycling  profile  was  repeated  continuously 
until  the  BSF-scaled,  3.4  kWh  of  energy  had  been  removed.  The 
cell  was  then  recharged  to  ~90%  SOC  and  allowed  to  rest  for  the 
times  indicated  in  Table  1.  The  discharge-rest-recharge  process 
was  repeated  continuously  for  32  days  at  the  test  temperature,  at 
which  time  an  HPPC  test  was  performed  to  re-characterize  the  cell 
and  to  gauge  capacity  fade  and  resistance  rise. 

The  data  from  the  HPPC  test  were  used  to  calculate  resistance 
vs.  voltage  curves.  The  discharge  and  regen  resistances  ( Rd  and  Rr, 
respectively)  were  calculated  using  the  equations  shown  in  Eq.  (1 ). 
For  Rd,  Vti  and  /ti  (see  Fig.  2  for  definitions  of  t\  and  t2)  are  the 
voltage  and  current,  respectively,  right  before  the  discharge  pulse, 
and  Vt2  and  Jt2  are  the  voltage  and  current  10  s  later.  For  Rr ,  Vt3,  Jt3, 
Vt4  and  /t4  (see  Fig.  2  for  definitions  of  t3  and  t4)  follow  analogously. 


Vtl  ~  Vt2  . 


kx  -  vt2  ’ 


Rr  = 


vt3  -  vt4 

U3  -  Vt4 


(1) 


An  example  plot  is  shown  in  Fig.  5  for  a  typical  cell.  As  the  cell 
aged,  the  discharge  and  regen  resistances,  generally,  increased  with 
time. 
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Fig.  5.  Discharge  and  regen  resistance  vs.  open  circuit  voltage  with  time  for  a  typical 
cell. 


Fig.  6.  The  change  in  relative  cell  capacity  vs.  time  for  the  calendar  life  cells.  The  val¬ 
ues  in  the  legend  are  aging  temperature  in°C;  the  negative  values  indicate  capacity 
loss.  The  cells  at  30,  40,  and  50  °C  tended  to  follow  a  different  capacity  loss  mech¬ 
anism  than  those  at  60  °C.  The  data  are  shown  as  markers,  and  the  fit,  as  a  solid 
curve. 


2.2.  Data  reduction 

The  cell  capacity  and  resistance  at  60%  SOC  were  collected  and 
normalized  to  t=0  before  curve  fitting.  Normalizing  the  data  from 
each  cell  to  that  point  minimized  the  effects  of  cell-to-cell  variation 
and  emphasized  the  relative  change  in  the  data  from  each  cell. 

For  most  of  the  work  described  below,  the  change  in  the  relative 
data  from  each  test  condition  was  fit  using  the  LINEST  function  in 
Microsoft®  EXCEL  or  using  the  non-linear  regression  function  in 
Systat®  SigmaPlot®.  In  some  cases,  such  as  in  CS  cycling  experiment 
that  used  SOC  as  the  experimental  parameter  (see  Section  3.2.2),  the 
change  in  the  normalized  data  was  fit  using  Battery  Life  Estimation 
(BLE)  software  [22]. 

Briefly,  the  BLE  software  is  based  on  the  mathematics  and 
reasoning  given  in  Ref.  [23a].  The  BLE  software  fits  a  proposed 
degradation  model  to  experimental  data  using  robust  linear  and 
non-linear  regression  techniques.  The  software  incorporates  a  sim¬ 
ple  error  model  that  attempts  to  quantify  the  level  of  deviations 
between  the  data  and  the  degradation  model.  It  also  uses  the  fitted 
model  to  estimate  battery  life  (Monte  Carlo  techniques).  The  non¬ 
linear  regression  uses  the  Levenberg-Marquardt  technique  [23b]. 
In  each  type  of  regression,  each  point  is  analyzed  and  iteratively 
re-weighted  using  the  biweight  transformation  [23c].  The  biweight 
transformation  is  most  commonly  applied  to  the  residuals  from  a 
linear  or  non-linear  fit.  The  biweight  transformation  has  the  advan¬ 
tage  that  it  allows  the  analyst  to  carry  out  subsequent  weighted 
linear  or  non-linear  fits  that  are  robust  and  resistive  to  outliers  in 
the  data  [23c]. 

3.  Results 

3.1.  Calendar  life  cells 

3.1.1.  Effect  of  temperature  (60%  SOC) 

Capacity.  The  change  in  the  relative  capacity  of  the  cells  versus 
time  is  shown  in  Fig.  6.  From  the  figure,  there  were  two  processes 
occurring.  One  capacity  fade  mechanism,  in  the  temperature  range 
of  30-50  °C,  followed  linear-with-time  kinetics,  and  the  other,  at 
60 °C,  followed  parabolic  (t1/2)  kinetics. 

From  the  fit  using  A(Rel.  capacity )  =  -exp(A  +  B/T)t,  where  A  and 
B  are  constants,  t  is  calendar  time,  and  T  is  absolute  temperature, 
the  data  from  the  30,  40,  and  50  °C  tests  also  followed  Arrhenius- 
like  kinetics.  The  fitting  parameters  and  the  standard  errors 
(s.e.)  for  this  group  of  data  were  A  =  3.84  x  1 01  (s.e.  =  4.86  x  1 0-1 ), 
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Fig.  7.  Change  in  relative  resistance  vs.  time  for  the  calendar  life  cells.  The  data  are 
shown  as  markers,  and  the  fit,  as  a  solid  curve. 

B=  —1.269  x  104  (1.55  x  102)I<,  and  r2  =  0.97.  From  the  high  value 
of  r2,  the  fit  was  very  good.  The  value  of  Ea,  the  activation  energy 
calculated  from  the  value  of  B,  for  this  degradation  process  was 
1 05.50  kj  mol-1. 

The  data  from  the  60  °C  were  fit  to  the  equation,  A(Rel.  capac- 
ity)  =  at1l2  +  c,  where  a  and  c  are  fitting  constants.  From  the  fit, 
the  values  of  a,  c,  and  r2  were  -8.33  x  10-1  (8.65  x  10-3)y-1/2, 
-1.44  x  10-1  (4.54  x  10-3),  and  0.99,  respectively.  Again,  the  high 
value  of  r2  indicated  that  the  fit  was  very  good.  The  nonzero  value 
of  c  reflects  the  error  in  the  average  capacity  measurement. 

Resistance.  The  values  of  cell  resistance  at  60%  SOC  were  treated 
as  described  above.  A  plot  of  the  change  in  the  relative  values  is 
given  in  Fig.  7.  From  the  plot,  all  data  sets  seemed  to  follow  linear- 
with-time  kinetic  laws. 

The  values  of  A  and  B  from  the  least-squares  fitting  were 
3.54  xlO1  (2.97  xlO-1)  and  -1.26  x  104  (^xlO^K  and  the 
value  of  r2  was  0.99,  indicating  that  the  fit  was  excellent.  The  esti¬ 
mated  value  of  Ea  was  96.10  kj  mol-1  for  this  degradation  process. 

3.1.2.  Effect  of  SOC  (30°  C) 

Capacity.  The  effect  of  cell  SOC  on  cell  capacity  is  shown  in  Fig.  8, 
along  with  the  least-squares  fit  to  the  data.  From  the  figure,  capacity 
decreases  with  time  at  all  SOCs.  Using  the  linear-with-time  kinetic 
expression,  At  +  B ,  the  data  were  fit  to  this  equation  with  values  of  r2 
of  0.99  or  better.  The  values  of  the  least-squares  fitting  parameters, 
along  the  standard  error,  for  these  curves  are  given  in  Table  2 


Fig.  8.  Change  in  relative  cell  capacity  vs.  time  and  SOC  for  the  calendar  life  cells. 
The  data  are  represented  as  markers  on  the  plot,  and  the  least-squares  fit,  as  solid 
curves.  The  numbers  in  the  legend  represent  the  SOCs  used  in  the  experiment. 


Table  2 

Fitting  parameters  from  the  calendar  life  SOC  cells. 


%soc 

A(s.e.),  y_1 

B(  s.e.) 

30 

-3.40  x  10“2  (9.62  x  10“4) 

1.82  x  lO-2  (4.74  x  10“4) 

60 

-4.59  xlO”2  (9.55  xlO”4) 

1.49  x  10-2  (4.71  x  10“4) 

90 

-8.74xl0-2  (2.32  xlO”3) 

1.32  x  lO"2  (1.14  xlO"4) 

From  the  data  in  Table  2,  the  values  A  seem  to  correlate  with 
SOC.  Indeed,  the  correlation  coefficient  between  the  values  of  SOC 
and  the  values  of  A  was  -0.95,  indicating  that  the  rate  of  capacity 
decline  increased  with  increasing  SOC. 

Resistance.  The  effects  of  SOC  and  time  on  the  resistance  of  the 
cells  are  shown  in  Fig.  9.  It  can  be  seen  from  the  figure  that  the 
effect  of  SOC  on  cell  resistance  was  not  simple.  At  30%  SOC,  the 
resistance  of  the  cells,  effectively,  did  not  change  during  the  course 
of  the  experiment;  on  the  other  hand,  the  resistances  of  the  cells 
at  60  and  90%  SOC  did  change.  In  the  60  and  90%  SOC  data,  the 
resistance  did  not  change  significantly  during  the  time  interval  of 
0.17  to  0.53y.  After  that  point  the  resistance  in  these  cells  tended 
to  increase. 

3.2.  CS  cycle  life 

3.2.1.  Effect  of  temperature  ( 60%  SOC) 

Capacity.  The  data  from  the  cycling-at-temperature  experiment 
are  shown  in  Fig.  10a,  along  with  fits  to  the  data.  The  data  could 
not  be  fit  to  an  Arrhenius-like  expression;  there  were  excessive 
errors  in  the  resulting  fit,  most  likely  from  the  fewer  number  of 
observations  at  50  °C.  Fewer  observations  at  the  high  temperature 
would  skew  the  fitting  because  of  the  uneven  weight  given  to  the 
observations  that  are  present  in  the  lower-temperature  studies. 
Consequently,  the  data  were  fit  as  four  independent  curves  to  the 
general  equation,  At?/2  +  B,  indicating  parabolic  kinetics.  The  values 
of  the  fitting  parameters  and  r2  are  given  in  Table  3 

An  estimate  of  the  activation  energy,  Ea,  for  this  degradation  pro¬ 
cess  was  calculated  from  the  values  of  A  in  Table  3.  From  the  slope 
of  the  line  of  ln(|A|)  vs.  1/T,  the  estimate  of  Ea  was  64.80  kj  mol-1 
(r2  =  0.99). 

Comparing  these  results  to  those  found  in  the  calendar  life  test 
shows  the  effects  that  CS  cycling  imposed  on  the  cells  (see  Fig.  1  Ob). 
In  all  cases,  CS  cycling  increased  capacity  fade  and,  in  most  cases, 
CS  cycling  changed  the  kinetic  behaviors.  Capacity  fade  in  the  cal¬ 
endar  life  test  tended  to  follow  Arrhenius-like,  linear-with-time 
kinetics  in  the  temperature  range  of  30-50 °C;  those  from  the  CS 


Fig.  9.  Change  in  relative  cell  resistance  vs.  time  and  SOC  for  the  calendar  life  cells. 
The  numbers  in  the  legend  represent  the  SOCs  used  in  the  experiment. 
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Table  3 

Values  of  the  fitting  parameters  and  r2  from  the  CS  cycling  at  temperature  tests. 
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Temperature,  °C 

A,  y-'112  (s.e.) 

B  (s.e.) 

r2 

30 

-8.23xl0-2  (2.55  xlO-3) 

2.37  x  10-2  (1.68  xl0“3) 

0.93 

40 

-2.10  x  10-1  (2.96  x  10-3) 

5.65  x  10-2  (1.94  x  10-3) 

0.99 

50 

-4.80  xlO-1  (2.63xl0-2) 

8.48x  10-2  (1.57x  10-2) 

0.94 

60 

-8.15  xlO-1  (1.33x10-2) 

9.07  x  IQ”2  (6.23  xl0“3) 

0.99 

cycle  test  exhibited  t1/2  kinetics.  Only  at  60  °C  in  the  calendar  test 
were  parabolic  kinetics  observed. 

Additionally,  cycling  increased  the  average  fade  rates.  At  40  °C, 
the  increase  was  modest  from  ~1 0%  to  ~1 2%  in  ~0.88y.  The  increase 
was  more  noticeable  at  60 °C,  from  35%  to  40%  in  ~0.35y.  Clearly, 
from  these  results,  the  majority  of  the  capacity  loss  was  due  to  the 
time  at  elevated  temperatures;  CS  cycling  had  a  relatively  minor 
influence  in  the  capacity  loss  process. 

Resistance.  The  changes  in  the  relative  resistances  from  the 
cycling-at-temperature  test  are  shown  in  Fig.  11,  along  with  fits 
to  these  data.  These  data  were  fit  to  the  equation,  A(Rel.  resis¬ 
tance )  =  exp(A  +  B/T)t ,  where  T is  absolute  temperature  and  t  is  time. 
The  values  for  A  and  B  were  33.72  (0.84)  and  -1 1047.04  K  (277.14), 
respectively.  From  the  value  of  r2,  0.99,  the  fit  was  very  good. 
The  calculated  curves  using  the  fitted  parameters  are  also  shown 
in  Fig.  11.  From  the  value  of  B ,  the  estimated  value  of  Ea  was 
91 .84  kj  mol-1  for  this  process. 

Comparing  these  fitting  results  with  those  obtained  from  the 
calendar  life  test  (see  Fig.  12)  shows  that  most  of  the  resistance 


Fig.  10.  (a)  Effect  of  temperature  of  CS  cycling  on  the  change  in  relative  capacity  vs. 
time.  The  numbers  in  the  legend  indicate  the  temperature  in°C.  Higher  tempera¬ 
tures  accelerated  capacity  loss,  (b)  Change  in  relative  capacity  vs.  time  showing  the 
difference  in  the  fitted  curves  from  the  calendar-  and  CS  cycle-life  tests.  CS  cycling 
increase  capacity  fade  and  can  change  the  kinetic  behavior  of  the  cells.  The  solid 
lines  represent  the  fitted  curves  from  the  calendar-life  test  and  the  dashed,  those 
from  CS  cycling. 
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Fig.  11.  Effect  of  temperature  on  the  change  in  relative  resistance  vs.  time.  The  num¬ 
bers  in  the  legend  indicate  the  temperature  in°C.  Higher  temperatures  accelerated 
resistance  increase.  The  observations  at  about  0.78y  from  the  40  °C  test  were  missing 
from  the  data. 

increases  are  due  to  calendar-life  effects  at  all  temperatures.  At  30 
and  40  °C,  there  were  no  significant  differences  between  the  cal¬ 
endar  and  CS  cycle  curves.  Only  at  50  and  60  °C  were  significant 
differences  observed;  the  cells  from  the  CS  cycle  test  displayed 
about  9-19%  additional  resistance  increase. 

3.2.2.  Effect  of  SOC  (30  °C) 

Capacity.  The  changes  in  relative  capacity  vs.  time  data  for  cells 
cycled  at  30,  60,  and  90%  SOC  are  given  in  Fig.  13,  along  with  the 
least-squares  fit  to  the  data.  Here,  the  change  in  relative  capacity 
was  fit  to  the  equation,  a  +  b  x  SOC  x  t,  where  a  and  b  are  fitting 
constants,  SOC  is  the  state  of  charge  expressed  as  a  decimal  frac¬ 
tion  and  t  is  time.  From  the  robust  regression,  the  values  of  a 
and  b  were  1.49  x  10"2  (8.01  x  lO”6)  and  -9.83  xl0“2  SOC-1  y-\ 
respectively.  Due  to  the  high  cell-to-cell  variability  (~11%),  the 
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Fig.  12.  Change  in  relative  resistance  vs.  time  comparing  the  fitting  results  from  the 
calendar  and  CS-cycling  tests.  The  numbers  next  to  each  set  of  lines  represent  the 
test  temperature.  The  lines  from  the  calendar  test  are  shown  in  solid  lines  and  those 
from  the  CS  cycle  test,  as  dashed  lines. 


10218 


J.  Belt  et  al /Journal  of  Power  Sources  196(2011)  10213-10221 


Time,  y 

Fig.  13.  Effect  of  SOC  on  the  change  in  relative  capacity  vs.  time.  The  numbers  in  the 
legend  indicate  the  SOC.  The  markers  represent  the  data  points  and  the  solid  lines, 
the  fitted  curves.  Higher  SOCs  tended  to  accelerate  capacity  loss. 


value  of  r2  was  low,  0.86.  The  results  of  this  experiment  are  con¬ 
sistent  with  those  described  in  the  calendar  life  tests  using  SOC  as 
the  experimental  variable,  that  rate  of  capacity  fade  increased  with 
increasing  SOC. 

Resistance.  The  effect  of  SOC  and  cycling  on  the  relative  resis¬ 
tance  of  the  cells  is  shown  in  Fig.  14,  along  with  fits  to  the  data.  As 
shown  in  Fig.  14,  the  change  in  relative  cell  resistance  was  complex. 
At  60  and  90%  SOC,  during  the  time  interval  of  0  to  ~0.44y,  the  resis¬ 
tance  of  the  cell  decreased.  After  ~0.5y,  it  increased.  However,  the 
data  from  the  30%  SOC  cells  shows  that  resistance  was  still  decreas¬ 
ing  during  the  time  of  the  test.  These  observations  imply  that  there 
was  a  voltage-driven  point  where  the  aging  process  fundamentally 
changed. 

It  is  interesting  to  note  that,  once  the  cell  was  operated  above  a 
certain  voltage,  the  inflection  point  in  the  resistance  vs.  time  data 
did  not  change.  That  is,  the  inflection  point  from  the  60  to  90%  SOC 
cells  was  at  the  same  point  in  time,  ~0.44y.  Beyond  the  inflection 
point,  the  trends  in  the  data  appear  to  be  the  same.  That  is,  there 
was  little  difference  in  the  trend  in  the  data  at  60  and  90%  SOC. 

3.3.  CD  cycle  life 
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Fig.  15.  (a)  Change  in  rel.  capacity  vs.  time  for  the  CD  cycle  life  test.  The  data  are 
indicated  as  markers,  and  the  fits,  as  solid  curves.  The  30-50  °C  cells  follow  one 
mechanism,  and  those  at  60  °C,  another,  (b)  Rates  of  change  for  two  kinetic  expres¬ 
sions,  x  =  t1/2  andx  =  lnt.  Initially,  d(lnt)/dt>d(t1/2)/dt. 


3.3 A.  Effect  of  temperature 

Capacity.  The  change  in  relative  capacity  from  this  experiment 
is  given  in  Fig.  15a.  From  the  figure,  the  effects  of  two  capacity¬ 
consuming  processes  are  seen.  Indeed,  comparing  Figs.  6  and  15a 
shows  that  they  appear  to  be  similar.  Both  display  the  effect  of 


Fig.  14.  Effect  of  SOC  on  the  change  in  relative  resistance  at  60%  SOC  vs.  time.  The 
numbers  in  the  legend  indicate  the  SOC.  All  cells  exhibited  an  interval  where  resis¬ 
tance  decreased.  This  interval  was  long  at  low  SOC  and  relatively  short  at  higher  SOC 
values. 


linear-with-time  processes  in  the  temperature  range  of  30-50  °C 
and  a  nonlinear  process  at  60  °C.  As  described  above  in  Section 
3.2.1,  the  data  could  not  be  fit  to  an  Arrhenius-like  expression; 
there  were  excessive  errors  in  the  resulting  fit,  most  likely  from 
the  fewer  number  of  observations  at  50  °C.  Fewer  observations  at 
the  high  temperature  would  skew  the  fitting  because  of  the  uneven 
weight  given  to  the  observations  that  are  present  at  longer  times. 
Fitting  the  30-50  °C  data  to  A(Rel.  capacity )  =  a t+b  yielded  the  val¬ 
ues  for  a,  b  and  r2  shown  in  Table  4.  Based  on  the  high  values  of  r2, 
these  fits  were  very  good 

At  first  look,  fitting  the  60  °C  data  indicated  that  equations  simi¬ 
lar  those  used  by  Broussely  et  al.  [24]  were  directly  applicable.  That 
is,  dx/dt  was  proportional  to  1/x,  where  x  is  the  concentration  of 
lithium  (i.e.,  cell  capacity).  Fitting  the  data  using  the  general  equa¬ 
tion,  t= k\X2  +  k2x,  yielded  a  reasonable  value  of  r2  of  0.93.  However, 
the  high  initial  rate  of  change  shown  in  the  data  was  not  truly  repre¬ 
sented.  The  data  were  better  described  by  using  a  mechanistically 
related,  lnt  expression.  Both  t1/2  and  lnt  can  be  used  to  describe 
diffusion-limited  processes  [25,26],  but  those  that  depend  on  lnt 
tend  to  change  faster  (dx/dt)  initially  (see  Fig.  15b).  Fitting  the  data 
to  the  expression  given  in  Eq.  (3)  yielded  the  following  parameters: 
D  =  -2.25  x  10-1  (4.77  x  1 0-3)  and  F  = -7.40  xlO-1  (7.54  xlO”3). 
The  value  of  r2  now  increased  to  0.99. 

A(Rel.  capacity)  =  D  In  t  +  E  (3) 
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Table  4 

Fitting  results  from  the  effect  of  temperature  during  CD  cycling  on  relative  capacity  test. 
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Temperature,  °C 

a,  y_1  (s.e.) 

b  (s.e.) 

r2 

30 

-2.28  xlO”1  (1.81  xlO-3) 

-1.49xl0-2  (9.84  xl0“4) 

0.99 

40 

-3.52  xlO-1  (3.08  xlO-3) 

-2.15  x  lO-2  (1.67  x  10-3) 

0.99 

50 

-8.30  xlO”1  (9.54  xlO-3) 

-3.05x  10-2  (3.26  xlO-3) 

0.99 

The  curves  in  Figs.  6  and  1 5a  show  the  net  effect  of  CD  cycling  on 
the  cells.  The  addition  of  cycling  increased  the  total  capacity  fade 
of  the  cells  as  well  as  the  average  rate  of  fade.  For  example,  adding 
cycling  increases  the  total  degradation  by  about  a  factor  of  2  in 
~0.53y  at  50  °C,  from  about  22%  to  about  45%.  The  effect  of  cycling 
on  the  60  °C  cells  can  also  be  seen,  but  the  net  effect  is  not  as  great. 
Here,  the  increased  loss  was  from  40%  to  about  55%  in  0.43y,  along 
with  a  change  in  the  fade  mechanism. 

Resistance.  The  resistance  data  from  the  CD  cycling  at  tempera¬ 
ture  tests  and  the  fits  to  these  data  are  given  in  Fig.  1 6.  An  inspection 
of  the  data  in  Fig.  16  shows  that  three  different  mechanisms  are 
present.  At  30  °C,  the  rate  of  change  increased  with  time;  at  40  and 
50 °C,  it  remained  constant;  and  at  60  °C,  it  tended  to  slow  with 
time. 

No  simple  kinetic  expression  was  found  to  fit  the  data  shown  in 
Fig.  16;  these  data  were  fit  independently  to  three  equations.  The 
data  from  the  30  °C  cells  was  fit  by  the  equation  of  the  form,  Atz  +  C. 
The  fitting  parameters  for  this  equation  were  0.71  (0.02 )y~  3 12, 1.52 
(0.10),  3.80  x  10-3  (1.19  x  10-3),  and  0.99  for  A,z,  C,  and  r2,  respec- 
tively.  Even  though  there  were  only  two  data  sets,  an  Arrhenius 
expression  (see  capacity  fade  in  this  section  for  form  of  the  equa¬ 
tion)  was  used  to  fit  the  data  from  the  40  and  50  °C  cells.  The  fitting 
parameters  here  were  12.00  (0.69),  3829.08  (219.30),  and  0.99  for 
A,  B,  and  r2,  respectively.  From  the  value  of  B,  the  estimated  value  of 
£a  was  3 1 .8  kj  mol-1 .  Finally,  the  data  from  the  60  °C  cell  test  was  fit 
using  a  parabolic  rate  law,  At1/2  +  C.  The  values  of  the  fitting  param¬ 
eters  were  1.13  (0.06)y_1/2,  -0.12  (0.03),  and  0.99  for  A,  C,  and  r2, 
respectively. 

Comparing  Figs.  7  and  1 6  shows  the  effect  of  CD  cycling  on  resis¬ 
tance  increase.  At  30 °C,  the  mechanism  changes  from  linear  to 
non-linear  and  the  average  resistance  changes  from  an  increase  of 
~10%  to  ~50%  in  0.79y.  It  is  interesting  to  note  that,  at  60  °C,  along 
with  the  mechanism  change,  the  average  rate  of  resistance  increase 
actually  decreased  from  ~110%  to  ~64%  in  0.43y. 

3.3.2.  Effect  of  rest  time  (30  °C) 

Capacity.  The  change  in  relative  capacity  of  the  cell  with  elapsed 
time  and  rest  time  is  shown,  along  with  a  fit  to  the  data,  in  Fig.  17. 
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Fig.  16.  Change  in  relative  resistance  vs.  time.  The  numbers  in  the  legend  represent 
the  test  temperatures.  The  data  are  indicated  as  markers,  and  the  fits,  as  curves. 
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Fig.  17.  Change  in  relative  capacity  vs.  time.  The  values  in  the  legend  indicate  the 
rest  time  between  charge  and  discharge  half-cycles  in  minutes.  The  data  are  shown 
as  markers,  and  the  fit  using  the  best  set  of  parameters,  as  the  solid  curves. 


Fig.  1 7  shows  that,  with  no  rest  time,  the  degradation  rate  was  the 
lowest.  With  longer  rest  times,  the  degradation  rate  increased.  The 
data  were  fit  using  AQ=  Ctz  +  Dtp,  where  AQ.  is  the  change  in  the 
relative  parameter  of  interest,  such  as  capacity  or  cell  resistance, 
t  is  calendar  time,  p  is  rest  time  expressed  in  calendar  years,  and 
C,  z,  and  D  are  fitting  constants.  The  fitting  constants,  C,  z,  and  D, 
were  -0.13  (1.80  x  10-3),  0.71  (1.87  x  \0~2),  and  -941.42  (51.74), 
respectively.  The  high  value  of  r2,  0.99,  indicated  that  the  fit  was 
very  good. 

The  effect  of  CD  cycling  on  capacity  decline  can  be  seen  by  com¬ 
paring  the  fitted  curves  from  the  calendar  life  test  and  Fig.  17;  this 
comparison  is  shown  in  Fig.  18.  As  expected,  capacity  loss  and  the 
average  rate  of  capacity  loss  were  greater  in  the  cycling  test  than 
those  from  the  calendar  test.  From  Fig.  18,  the  capacity  loss  in  the 
O-min-rest-time  CD  cycling  test  was  about  9.3%  in  ~0.8y,  while  that 
from  in  the  30  °C  calendar  test  was  about  3%  in  the  same  amount 
of  time.  The  3-fold  increase  in  capacity  fade  reflects  the  additional 
stress  caused  by  CD  cycling. 

The  relative  magnitudes  of  C  and  D  indicate  the  influences  that 
elapsed  time  and  rest  time  have  on  capacity  fade.  In  the  cycling 


Fig.  18.  Fitted  curves  from  Fig.  17  and  the  calendar  data  for  cells  tested  at  30  °C, 
showing  the  comparison  of  calendar  aging  with  that  from  CD  cycling.  As  expected, 
CD  cycling  greatly  increased  the  rate  of  capacity  loss. 
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Fig.  19.  Change  in  relative  resistance  vs.  calendar  time.  The  values  in  the  legend 
reflect  the  rest  time  in  minutes.  The  data  are  shown  as  markers,  and  the  fit,  as  a 
solid  curve. 

experiment,  it  appears  that  there  is  a  strong  interaction  between 
the  total  amount  of  rest  time  between  charge  and  discharge  and 
the  duration  of  the  test.  As  shown  in  Fig.  9,  capacity  fade  is  strongly 
influenced  by  time  at  high  SOC.  With  increasing  rest  time,  the  cells 
spend  more  time  idle  at  high  voltage  (90%  SOC),  where  degradation 
would  be  expected  to  be  faster. 

In  the  calendar  test  discussed  above,  capacity  fade  at  30  °C  was 
proportional  to  t.  With  CD  cycling,  the  mechanism  became  propor¬ 
tional  to  t0-71 . 

Resistance.  The  change  in  the  relative  cell  resistance,  along 
with  fits  of  these  data,  is  shown  in  Fig.  19.  From  the  figure,  the 
change  in  cell  resistance  increases  geometrically  with  elapsed  time 
and  is  directly  proportional  to  rest  time.  The  data  were  fit  using 
A R  =  Ctz  +  Dtp  +  E ,  where  A R  is  the  change  in  the  relative  resistance, 
t  is  calendar  time,  p  is  rest  time  expressed  in  calendar  years,  and  C, 
z,  D  and  E  are  fitting  constants 

The  fit  of  these  data  using  the  same  equation  used  in  the 
capacity  fade  discussion  in  Section  3.3.2  further  amplifies  these 
observations;  the  values  of  the  fitting  parameters,  C,  z,  D,  E 
and  r2  are  1.43  x  10-1  (3.10  x  10-3),  1.62  (3.30  x  10-2),  1.97  x  104 
(3.32  x  101),  and  0.99,  respectively. 

Fig.  20  compares  this  fit  to  that  from  the  calendar  life  data  at 
30  °C.  Again,  cycling  increased  the  change  in  cell  resistance  as  well 
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Fig.  20.  Comparison  of  the  fitted  curves  from  the  rest-time  tests  and  that  from  the 
30  °C  calendar  life  test.  The  curve  representing  the  30  °C  calendar  life  fit  has  been 
artificially  translated  0.073  units  down  to  make  the  comparison  between  the  curves 
easier. 


as  the  average  rate  of  change.  For  example,  the  resistance  of  the 
cells  with  no  rest  time  increased  1 0.5%  in  ~0.8y;  while,  in  the  30  °C 
calendar  life  test,  it  increased  only  about  5.7%  in  the  same  amount 
of  time.  The  ~2-fold  increase  in  resistance  reflects  the  additional 
changes  caused  by  CD  cycling. 

As  discussed  above,  the  relative  magnitudes  of  the  fitting  con¬ 
stants,  C  and  D,  indicate  the  influence  that  calendar  time  and 
rest  time  have  on  resistance  change.  In  the  cycling  experiment,  it 
appears  that  there  is  a  strong  interaction  between  the  total  amount 
of  rest  time  between  charge  and  discharge  and  the  duration  of  the 
test. 

4.  Discussion 

The  capacity  fade  and  resistance  increase  results  from  the  cal¬ 
endar  life  test  in  the  temperature  range  of  30-50  °C  and  60%  SOC 
appeared  to  follow  Arrhenius-like  kinetics.  Similarly,  resistance 
increase  seemed  to  follow  Arrhenius-like  kinetics.  Thus,  capacity 
fade  and  resistance  increase  may  be  caused  by  thermally-activated 
processes  which  were  linear  with  time. 

It  was  interesting  to  note  that  resistance  increase  followed  two 
processes  in  this  test,  one  that  was  active  in  the  temperature  range 
of  30-50  °C,  and  the  other,  at  60  °C.  The  rate  of  resistance  increase 
at  60 °C  did  not  depend  on  t;  instead  it  depended  on  t1/2.  At  this 
temperature,  resistance  increase  may  be  due  to  the  growth  of  a 
solid  electrolyte  interface  layer  (SEI),  the  growth  of  which  can  be 
proportional  to  t1/2  [13,27,28]. 

Increasing  or  decreasing  the  SOC  also  had  effects  on  the  rates 
of  capacity  fade  and  resistance  increase.  The  rate  of  capacity  fade 
tended  to  increase  with  increasing  SOC.  That  is,  time  at  high  volt¬ 
age  was  deleterious  to  capacity  retention.  On  the  other  hand,  the 
changes  in  resistance  were  not  as  simple;  cell  resistance,  effectively, 
did  not  change  at  30%  SOC.  At  60%  and  90%  SOC,  cell  resistance 
did  increase,  but  there  was  no  significant  difference  in  its  behavior 
with  time  at  these  SOCs.  This  observation  implies  that,  above  30% 
SOC,  the  rate  of  resistance  increase  was  not  sensitive  to  SOC  (or  cell 
voltage). 

The  results  from  CS  cycling  showed  that  the  rates  of  perfor¬ 
mance  change  were  very  similar  to  those  observed  in  the  calendar 
life  experiment.  On  the  other  hand,  the  results  from  the  CD  cycling 
test  were  not.  There  were  dramatic  changes  in  the  mechanism(s) 
and  rates  of  performance  change  in  the  CD  cycling  test.  Since 
the  amount  of  energy  removed  during  CD  cycling  differs  from 
that  during  CS  cycling  by  about  two  orders  of  magnitude,  (from 
Figs.  3  and  4, 3.8  and  545.61  Wh  per  profile,  respectively),  the  rates 
of  performance  change  shown  above  in  CS  and  CD  cycling  may  be 
linked  to  energy  removal  or  the  wider  voltage  range  used  during 
CD  cycling. 

5.  Conclusions 

One  hundred-seven  commercially-available,  off-the-shelf,  1.2- 
Ah  cells  were  tested  for  calendar  life  and  CS  cycle-  and  CD  cycle-life 
using  the  new  USABC  PFIEV  Battery  Test  Manual.  Here,  the  effects  of 
temperature  on  calendar  life,  on  CS  cycle  life,  and  on  CD  cycle  life; 
the  effects  of  SOC  on  calendar  life  and  on  CS  cycle  life ;  and  the  effects 
of  rest  time  on  CD  cycle  life  were  investigated.  The  results  indicated 
that  the  test  procedures  caused  performance  decline  in  the  cells 
in  an  expected  manner,  calendar  <  CS  cycling  <  CD  cycling.  In  some 
cases,  the  kinetic  law  changed  with  test  type,  from  linear-with-time 
to  about  t2.  Additionally,  temperature  was  found  to  stress  the  cells 
more  than  SOC,  causing  increased  changes  in  performance  with 
increasing  temperature. 

For  maximum  cell  life,  the  rest  time  during  CD  cycling  should 
be  minimal.  This  may  directly  impact  the  design  of  test  profiles. 
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